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k-cKP hierarchy, (t^, Tb) and (7^, as) matrix hierarchies. (2+l)-BDk-cKPH con- 
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. Dressing Method is also proposed for construction of exact solutions for equations 

from (2+l)-BDk-cKPH. As an example the exact form of multi-soliton solutions for 
^ ■ vector generalization of the Davey-Stewartson system is given. 

^ ■ PACS numbers: 05.45 Yv, 05.45-a, 02.30 Jr, 02.10 Ud 

cn 
o 



''^ Electronic mail: alex.chvartatskyy@gmail.com 
^^Electronic mail: y_sydorenko@franko.lviv.ua 



I. INTRODUCTION 



In the modern theory of nonhnear integrable systems, algebraic methods play an im- 
portant role (see the survey in^). Among them there are the Zakharov-Shabat dressing 
method^"'*, Marchenko's method^, and an approach based on the Darboux-Crum-Matveev 
transformations^'^. Algebraic methods allow us to omit analytical difficulties that arise in the 
investigation of corresponding direct and inverse scattering problems for nonlinear equations. 
A significant contribution to such methods has also been made by the Kioto group In 
particular, they investigated scalar and matrix hierarchies for nonlinear integrable systems 
of Kadomtsev-Petviashvili type (KP hierarchy). 

The KP hierarchy is of fundamental importance in the theory of integrable systems and 
shows up in various ways in mathematical physics. Several extensions and generalizations 
of it have been obtained. For example, the multi-component KP hierarchy contains several 
physically relevant nonlinear integrable systems, including the Davey-Stewartson equation, 
the two-dimensional Toda lattice and the three-wave resonant interaction system. There are 
several equivalent formulations of this hierarchy: matrix pseudo-differential operator (Sato) 
formulation, r-function approach via matrix Hirota bilinear identities, multi-component free 
fermion formulation. Another kind of generalization is the so-called "KP equation with self- 
consistent sources" (KPSCS), discovered by Melnikov^"— . In^^"— , k-symmetry constraints 
of the KP hierarchy which have connections with KPSCS were investigated. The resulting 
k-constrained KP (k-cKP) hierarchy contains physically relevant systems like the nonlinear 
Schrodinger equation, the Yajima-Oikawa system, a generalization of the Boussinesq equa- 
tion, and the Melnikov system. Multi-component generalizations of the k-cKP hierarchy 
were introduced in^. In the papers^i^^ the differential type of the gauge transformation 
operator was applied to the k-constrained KP hierarchy at first. In— differential and inte- 
gral type of the gauge transformation operators were applied to k-cKP hierarchy. In^^ the 
binary Darboux transformations for dressing of the multi-component generalizations of the 
k-cKP hierarchy were investigated. A modified k-constrained KP (k-cmKP) hierarchy was 
proposed irt^ii^i^. It contains, for example, the vector Chen-Lee-Liu and the modified KdV 
(mKdV) equation. Multi-component versions of the Kundu-Eckhaus and Gerdjikov-Ivanov 
equations were also obtained in^-, via gauge transformations of the k-cKP, respectively the 
k-cmKP hierarchy. 
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Moreover, iiM-^, (2+l)-dimensional extensions of the k-cKP hierarchy were introduced 
and dressing methods via differential transformations were investigated. IdM^, exact so- 
lutions for some representatives of the (2+l)-dimensional k-cKP hierarchy were obtained 
by dressing binary Darboux transformations. This hierarchy was also rediscovered recently 
in2^. Dressing methods via differential transformations for this hierarchy and its modified 
version were investigated in'^^. The (2+l)-dimensional k-cKP hierarchy in particular con- 
tains the DS-III (A; = 1), Yajima-Oikawa {k = 2) and Melnikov {k = 3) hierarchies. The 
corresponding Lax representations of this hierarchy (see ([9])) consists of one differential and 
one integro-differential operator. Our aim was to generalize the (2+l)-dimensional k-cKP 
hierarchy to the case of two integro-differential operators in Lax pair. It is essential 
to call this new hierarchy a bidirectional generalization of (2+l)-dimensional k-cKP hierar- 
chy or simply (2+l)-BDk-cKP hierarchy (see (l30l) ). We will consider this hierarchy in the 
most general matrix case. Another aim of our investigation is a construction of the Binary 
Darboux Transformation Dressing Method for (2+l)-BDk-cKP hierarchy ( l30l) . 

This work is organized as follows. In Section 2 we present a short survey of results 
on constraints for KP hierarchies and their (2+l)-dimensional generalizations. In Section 
mil we introduce (2+l)-BDk-cKP hierarchy. Members of the obtained hierarchy are also 
listed there. (2+l)-BDk-cKP hierarchy contains matrix generalizations of Davey-Stewartson 
hierarchy (first members of it are two different matrix versions of DS-III system: DS-III- 
a ( HHj) and DS-III-b (149|) ). new Yajima-Oikawa and Melnikov hierarchies. In Section IIVI 
we consider dressing via binary Darboux transformation for the hierarchy constructed in 
Section IIIII Exact forms of solutions for some members of this hierarchy are also presented 
here in terms of Grammians. In the final section, we discuss the obtained results and 
mention problems for further investigations. Some examples of Lax representations from 
the bidirectional generalization of (2+l)-dimensional k-constrained modified KP hierarchy 
((2+l)-BDk-cmKP hierarchy) are also presented there. 



II. fc-CONSTRAINED KP HIERARCHY AND ITS EXTENSIONS 

To make this paper self-contained, we briefly introduce the KP hierarchy^, its k-symmetry 
constraints (k-cKP hierarchy), and the extension of the k-cKP hierarchy to the (2+1)- 
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dimensional case^*^. A Lax representation of the KP hierarchy is given by 

Li„ = [S„,L], n>l, (1) 

where L = D + UiD^^ + U2D^'^ + . . . is a scalar pseudodifferential operator, ti := x, D := 
and Bn := {L"-)+ '■= (-^")>o = -D" + J2i=o UiD^ is the differential operator part of L". The 
consistency condition (zero-curvature equations), arising from the commutativity of flows 

Bn,tk — Bk,tn + [Bn, Bk] = 0. (2) 

Let denote the formal transpose of Bn, i.e. BJ^ := {—l)"'D"' + '^^^Q{—iyD^uJ, where 
^ denotes the matrix transpose. We will use curly brackets to denote the action of an 
operator on a function whereas, for example, i?„ q means the composition of the operator 
Bn and the operator of multiplication by the function q. The following formula holds for 
Bnq and Bn{q}- Bn{q} = BnQ — (-B„g)>o. The k-cKP hierarchy^^"— is given by 

Lt^ = [Bn,L], (3) 

with the k-symmetry reduction 

Lfe := = Bk + qD-h. (4) 

The hierarchy given by ©-(jl]) admits the Lax representation (here G N is fixed): 

[Lk, Mn] =0, Lk = B, + qD-\ M„ = dt„ - Bn- (5) 

Lax equation ([5]) is equivalent to the following system: 

[L,, Mn]>o = 0, Mn{q} = 0, M;{r} = 0. (6) 

Below we will also use the formal adjoint B* := B^ = (-1)"D" + Er=o^(-l)'^'< of 
Bn, where * denotes the Hermitian conjugation (complex conjugation and transpose). In2^, 
multi-component (vector) generalizations of the k-cKP hierarchy were introduced, 

m m 

Lk:=L' = Bk + J2Yl l^^^jD-'rj = Bu + qMoD-'r'^ , (7) 

i=i j=i 

where q = {qi, . . . ,qm) and r = (ri,...,rm) are vector functions, A^o = ("^«j)i^=i is a 
constant m x m matrix. Aio can be annihilated by the change of functions: qA^o 
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j,T _^ j,T q _v. A^qI""'' — )■ r"""). However, we remain this matrix because it plays 
an important role in Section HVl where we consider dressing methods for (2+l)-BDk-cKP 
hierarchy (130|) . We shall note that extensions of hierarchies given by ([3]) and ([7]), namely 
(2+l)-dimensional k-cKP ([JU]) and (2+l)-BDk-cKP ([10]), (USD hierarchies are represented 
analogously to (E]) (with more general Lj. and M„ operators). 

For k = 1, the hierarchy given by ([3]) and ([7]) is a multi-component generalization of the 
AKNS hierarchy. For k = 2 and k = 3, one obtains vector generalizations of the Yajima- 
Oikawa and Melnikov^ii^ hierarchies, respectively. 

jj^^jSi^ a k-constrained modified KP (k-cmKP) hierarchy was introduced and investi- 
gated. Its Lax representation has the form 

[Lfc, M„] =0, Lk = Bk + qMoD^h'^D, M„ = dt„ - 5„. (8) 

where = + ^^Z^WjD^ ■ For k = 1,2,3, this leads to vector generalizations of the 
Chen-Lee-Liu, the modified mult i- component Yajima-Oikawa and Melnikov hierarchies. 

An essential extension of the k-cKP hierarchy is its (2+l)-dimensional generalization^i^, 
given by 

Lk = Pkdr, -Bk- qMoD-h^, M„ = arA^ - A^, 

Bk = D^ + Y!;zI UjD^, An = D^ + Y.Zo (9) 

where Uj and Vi are scalar functions, q and r are m-component vector-functions. Lax 
equation [L^, M„] = is equivalent to the system: 

[Lfc, M„]>o = 0, M„{q} = 0, Ml{v} = 0. (10) 

System ( fTOj) can be rewritten as 

dnBku = PkAn,r, + [^n, ^fc] + ([^n, qA^o^^^'r^] ) >o , ^n{q} = 0, M^jr} = 0. 
We list some members of this (2+l)-dimensional generalization of the k-cKP hierarchy: 
1. A; = 0, n = 2. 

Lo = f^odro - qMoD-'r^, M2 = a2dt, - - vq. (11) 
The commutator equation [Lq, M2] = is equivalent to the system: 

«2qt2 = Qxx + ^^oq, -o:2rJ^ = rj^ + r'^^^o, /3oVo,ro = -2(qAlor'^)x- (12) 
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After the reduction /3o G M, a!2 £ i^, r = A^o = -^O' ''^o = "^o, the operators Li 
and M2 in (fTT]) are skew-Hermitian and Hermitian, respectively, and (fT2|) becomes the 
DS-III system: 

"2qt2 = q^x + ^'oq, ^o^^o,ro = -2(qA^oq*)x- 

2. A; = 1, n = 2. Then has the form 

i^i = f3idr, -D- ciMoD-^^, M2 = a2dt, - - vo, (14) 
and the equation [Li, M2] = 0, is equivalent to the system, 

"2qt2 = (Ixx + voq, Q!2rt2 = -r^a; - VqT, 

[lb) 

(3iVo,Ti = Vo,x - 2(qA^or^)x- 

After the reduction /3i G M, a2 G iR, r = q, Mq = M^, Vq = vq, we obtain the DS-III 
system in the following form: 

a2qt2 = <ixx + Voq, /3iVo,ri = Vq^^ - 2(qA^oq*)x- (16) 

In Remark [1] we will show that DS-III systems f lT3|) and (fT6|) as well as all the other 
equations of the hierarchy related to the operators Lq and Li in ([9]) are equivalent 
via linear change of independent variables. 

3. A; = 1, n = 3. Now (j9]) becomes 

Li = ^,dr, -D- qMoD-h'^, M3 = asdt, - - v,D - ^o- (17) 

After the additional reduction q;3,/?i G M, ^Ao = A^q, Vi = Vi, Vq + Vq = Vix, the 
operators Li, M3 in (fT7|) are skew-Hermitian, and the Lax equation [Li,M3] = is 
equivalent to the following (2+l)-dimensional generalization of the mKdV system: 

tt3qt3 = ^xx + Viq^ + Voq, 
PiVo^Ti = fo,x - 3(qa;A^oq*)x, /3i^^i,ri = t'l.x - 3{qMoq*)x- 

This system admits the real version {Ai^ = M.q, q* = q""") 

ttsqta = (ixxx + viqx + lvi,xq, f^ivi^n = vx,x - 3(qAloq'^)x- (19) 
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4. A; = 2, n = 2. ([9]) takes the form 

L2 = Mr, - D'' - uo - (iMoD-\'^ , M2 = a2dt, - - uq. (20) 

Under the reduction a2, /32 G i^, r = q, -M-o = —-M.^, Wq = Uo '■= 2if, the operators L2 
and M2 in (12U|) become Hermitian and the Lax equation [L2, M2] = is equivalent to 
the following system: 

«2qt2 = flxx + 2Mq, 02^^*2 = /32^^r2 + (q-Moq*)a;, (21) 
which is a (2+l)-dimensional vector generalization of the Yajima-Oikawa system. 

5. k = 2, n = 3. Now IMl becomes 



L2 = P2dr, -D^- 2no - qTWo^-'r^, ^^2) 
M3 = a^dt, -D'- 3uoD - I {uo,^ + P2D-^{uo,r,} + qA^or^) • 

In formula f l22|) . D~^{iio,T2} denotes indefinite integral of the function UQ^r, with respect 
to X. With the additional reduction P2 G ^1^, «3 G M, lio = Uq := m, A^o = ^AIq 
r = q, this is equivalent to the following generalization of the higher Yajima-Oikawa 
system^ (in scalar case (m = 1) Manakov LAB-triad for this equation was proposed 
inM): 

«3qt3 = qxxx+ Swq^. + 1 (M:r + i32D-^{urJ + qA^oq*) q, ^23) 
[a^ut^ - ju^^^ - 3uu^ + f (qAloq^ - qxAloq% -|/32 (qAloq*)^J^ = f/32"r2r2. 

6. k = 3,n = 3. takes the form 

L3 = /33dr,-D^-uiD-^ui,,-qMoD-'r'^, M3 = asdt,- -uiD-^m,,. (24) 

With the additional reduction a^, G M, "Ui = Mi := u, Aio = AIq ~ Aio and 
r = q = q, this is equivalent to the following (2+l)-dimensional generalization of the 
Drinfeld-Sokolov-Wilson equation^"—: 

«3qt3 = ^ccxx + uqx + ^^ixq, 03^*3 = l^sUrs + 3(qAloq'^)x- (25) 

The following remark establishes connections between equations represented by Lax pairs 
(Lq, Mfc) and (Li, M^), {L^, M^) and (Li, M^) respectively (see formulae ([9])). 



Remark 1. The change of variables x := x — -^Tq, Ti '■= ^q'^o> •= q(3^)''o)? 

Vo{x,Ti) := Vo{x,tq) maps DS-III equation fT3\} to another form fT^) of it. The same 
change maps the "higher" equations related to operator Lq to the equations related to Li (see 
formulae 

Remark 2. The change of variables 

x:= x- j^Tk,ik ■= tk + f^n,ri := |^rfc, ffc_2(x, 4, n) := Uk-2{x,tk,Tk), 

~ _ ~ ^ 
q(x,tfc,ri) := q{x,tk,Tk), 

maps (2+1 )-dimensional generalization of Yajima-Oikawa system l[EJ\) to DS-III system / f7^) 
in the case k = 2. In the case k = 3 it maps (2+1)- extension of Drinfeld-Sokolov- Wilson 
equation [2B) to the real version of (2+1)- extension of mKdV equation [W) . An analogous 
relation holds between all other equations with Lax pairs (Lk, Mk) and (Li, Mk) where the 
operators Li, Lk and Mk are defined by formulae 

Thus, for A; = 1 we have the DS-III hierarchy (its first members are DS-III {k = 1, n = 2) 
and a special (2+l)-dimensional extension of mKdV (A; = l,n = 3), see (fT6|) and (fT8|) . 
For A; = 2, = 3, we have (2+l)-dimensional generalizations of the Yajima-Oikawa (in 
particular, it contains (1?!]) and (125]) ) and the Melnikov hierarchyi^, respectively. 



III. BIDIRECTIONAL GENERALIZATIONS OF (2+l)-DIMENSIONAL 
K-CONSTRAINED KP HIERARCHY 

In this section we introduce a new generalization of the (2+l)-dimensional k-constrained 
KP hierarchy given by ^ to the case of two integro-differential operators. One of them (the 
operator Mn^i in (130|) ) generalizes the corresponding operator M„ ([9]) and depends on two 
independent indices / and n. It leads to generalization of (2+l)-dimensional k-cKP hierarchy 
([9]) in additional direction / (/ = 1, 2, . . .). We do not consider the case I = because of the 
Remark [21 For further purposes we will use the following well-known formulae for integral 
operator hiD~^h2 constructed by matrix-valued functions hi and /i2 and the differential 
operator A with matrix-valued coefficients in the algebra of pseudodifferential operators: 

Ah{D'^h2 = {Ah{D-^h2)>o + A{hi]V-^h2, (27) 
h{D~^h2A = {hiV~'h2A)>o + hV~^[A^{h'^}]^, (28) 
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hiV-^hihsV-^hi = hiD-^{h2hi]V-^hi - hiV'^ {h2h:i}hi. (29) 

Consider the following bidirectional generalization of (2+l)-dimensional k-constrained 
KP hierarchy ([H]): 

Lk = Pkdrj, - Bk- qMoD'h^, Bk = Yl'j=o UjD^,Uj = Uj{x, r^, t„), G C, 

= a„at„ - A„ - Q E -=0 qb1A^oI5-'rT[; - j'], 1 = 1,... (30) 

= YJ'i=o'"iD\Vi = Vi{x,Tk,tn),an G C, 

where Uj and ^Jj, q and r are N x N and N x M matrix functions respectively; q[j] and r[j] 
are matrix functions of the following form: q[j] := {Lk)-'{q}, r^[i] := {{Ll)^ {r})~^ . 
The following theorem holds. 

Theorem 1. The Lax equation [Lk,Mn,i] = is equivalent to the system: 

[Lk, M„,,]>o = 0, M„,,{q} = Q(Lfc)'+i{q}, M;,{r} = Q(L^)'+i{r}. (31) 

Proof. From the equality [Lfc,M„_;] = [Lfc,M„j]>o + [Lfc,M„^;]<o we obtain that the Lax 
equation [Lk, Mn^i] = is equivalent to the following one: 

[Lk, Mn,i]>o = 0, [Lfc,M„,i]<o = 0. (32) 

Thus, it is sufficient to prove that [Lk, Mn,i]<o = <^==^ {Mn,i — Ci{Lky^^){q} = {M^i — 
Q(L^)'^^){r} = 0. Using bi-linearity of the commutator and explicit form (l30l) of the 
operators Lk and M„ ; we obtain: 



[Lk, M„,i]<o = Q E •=o[qb1-^oI^-'rT[/ _ ^-j^ ^^^^^ _ 5^]^^^ 

(33) 



l<0- 



After direct computations of each of the three items on the right-hand side of formula 
331) we obtain: 



1. 



ci EU^^[j]MoD-'r^[l - j],Pkdr, - BkU = -ci EU i^'^'lrAj] - BkMj]}) ■ 
■MoD-h^[l - j] - ci E5=o ^[j]MoD-' {^rl [I - j] + Bl{r^[l - j]}) . 

(34) 

Equality is a consequence of formulae (|27|) - (128|) . 
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2. 



ci T!,=MM,D-^r^ ^[j\M,D-h'^[l - j]]<o = ci EU ciMoD-'{r^q[j]} 



X 



X 



MoD-^r-^il - j] - Q E;=o qMoD-^D-^{v^c^[j]}Mor^[l - j]- 



(35) 



(36) 



Formula f l55]) follows from 

3. 

[ttn^t^ - An, qMo^'^rT]<o = (anqt„ - ^n{q}) A^oI^~^r^+ 

The latter equality is obtained via (l27!) - (l28|) . 
From formulae (l33|) - fl36|) we have 

[i^fc,M„,,]<o = Q (E5=oqb1-^o^-'^I{r[/ - j]} - E;=o^Uqb1}-^o^-'rT[; - + 

+M„,;{q}A^oI5-^rT - qA^oi?-^(M;,{r})T = q E5=o qb1-^oI5-^rT[/ - j + 1]- ^^^^ 
-Q Ej=o qb + AM^D-\^\l - 3\ + M„,^{q}A^oI^-'rT - qMo^^-^M^^jr})^ = 
= (M„, - Q(Lfc)'+i){q}A^oI?-'rT - qA1o^-i((M;, - c,{Uj^'){v}y . 

From the last equality we obtain the equivalence of the equation [L^, M„ = and (EI])- □ 
New hierarchy (l30l) consists of several special cases: 

1. = 0, Q = 0. Under this assumption we obtain Matrix k-constrained KP hierarchy^i. 
We shall also point out that the case = and q 7^ also reduces to Matrix 
k-constrained KP hierarchy. 

2. a„ = 0. In this case we obtain a new (l+l)-dimensional Matrix k-constrained KP 
hierarchy^'''. Their members are stationary with respect to tn members of (2+l)-BDk- 
cKP hierarchy ( l30i) . 

3. Q = 0, = 1, w„ = Ufc = 1, Vn-\ = Uk-i = 0. In this case we obtain (2+l)-dimensional 
k-cKP hierarchy given by ([9]). 

4. n = 0. In this case the differential part of Mq i in (130|) is equal to zero: Aq = 0. We 
obtain a new generalization of DS-III hierarchy. 
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5. Q = 0. We obtain (t^, rB)-Matrix KP Hierarchy that was investigated 

6. If we set / = we obtain (7^, o"B)-Matrix KP hierarchy that was investigated in^. 
However, we do not include / = in formula ( !30|) because of the following remark. 

Remark 3. In the case I = we obtain: 

Mn,0 = andt„ - An- CqCiMoD'^Y^ = {andt„ - CoPkdrJ - {An - CoBk) + CoLk (38) 

The last summand in the latter formula can he ignored because commutes with itself. 
Thus, we obtained an evolution differential operator again. It means that in the case 
I = the hierarchy given by the operators Mnfi and Lk after the change of independent 
variables (andt„ — CQ/3kdr^. — )■ (kndi^J coincides with hierarchy ^ in the matrix case. 

The following corollary follows from Theorem [1] 
Corollary 1. The Lax equation [Lfc,M„,^;] = 0, where 

Mn,l = Mn,i - CliLkY^' (39) 

and the operators and Mn are defined by / fS'Z^) . is equivalent to the system: 

[Lk, Mn,i]>o = 0, MnA^} = 0, ML{r} = 0. (40) 

Proof. It is evident that the operator commutes with the natural powers of itself: 
[Lfc, [LkY^^] = 0, / G N. Thus by bilinearity of the commutator we obtain that [L^, M„j] = 
holds if and only if [Lk,Mn,i\ =0. It remains to use Theorem [T] in order to complete the 
proof. □ 

As a result of the hierarchy given by operators (l30l) and a bilinearity of the commutator 
we obtain the following essential generalization of (130|) : 

Lk = (3kdrk - Bk- qA^o^""^r^, Bk = Yl'j=o UjD^,Uj = Uj{x, Tk, t„), f3k E C, 

Pn,m = " A„ - ZT=i ^1 {{eU cL[j]MoD-^r^ [I - j]) + (Lfc)'+i) , m = 1, 2, . . . 

An = YH=QViD\Vi = Vi{x,Tk,tn),an E C, 

(41) 

and the following corollary that is immediate consequence of Theorem [1) 
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Corollary 2. The commutator equation [Lk, Pn,m] = is equivalent to the following system: 

[Lk, Pn,m]>0 = 0, Pn,M = 0, PI^U} = 0. (42) 



For further convenience we will consider the Lax pairs consisting of the operators Lk ( l30|l 
and Mn^i ( !39|) (the operator M„ ^ is involved in equations for functions q and r; see formulae 
( 140|) ). Consider examples of equations given by operators Lk ( !30|) and Mn,i (139|) that can be 
obtained under certain choice of {k, n, I). 

I. k = 0. 

In this case the operator Lq fl3Up has the form: 

Lo = Podro - qMoD-'r'^. (43) 
For further simplicity we use the change of variables tq := y, (3 := /3o: 

Lo = I3dy - ciMoD-h'^. (44) 

1. n = 2,l = l. 

M2,i = M2,i - Ci(Lo)2 = a2dt, - c{D^ + v^) - c^^'^dl + 2ci/3q7WoI^-'rJ+ ^^^^ 
+2cil3ciMQD-^Y'^ dy, ceC. 



If c = we obtain that M2^i = Mq^i. We use c in formula (145|) in order not to consider 
separately the case n = 0. The commutator equation [Lo,M2^i] = is equivalent to 
the system: 

a2qt2 = cqxx + Ci/3'^qyy + cvoq + CiqA^o^', 
-a2rl = crj^ + Ci/J^rj^^ + cr^vo + ci^A^or^, (46) 
/3voy = -2(q7Wor^)x, = -2f3{r'^q)y. 
The Lax pair Lq f lH|) and M2^i fH5]) and the corresponding system (H^ were investi- 
gated in^"^. Consider additional reductions of pair of the operators Lq (jB]) and M2,i 
(li5|) and system (jl6]). After the reduction c, Ci G M, /3 G M, 0^2 G r""" = q*, 
A^o = -^0' operators Lq and M2,i are skew-Hermitian and Hermitian respectively, 
and (146|) takes the form 

a2qt2 = cqa;a; + Ci/3^qj^y + cvoq + CiqA^o^', ^^^^ 
f3voy = -2(qA^oq*)x, = -2(3{q*q)y. 

This has the following two interesting subcases: 
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(a) Ci = 0. Then we have 

a2qi2 = cq^x + cvoq, /3voy = -2(qA^oq*)x- (48) 

(b) c = 0. Then (H7|) takes the form 

a2qt, = Ci^^Hyy + CiqMoS, = -2/3(q*q)y. (49) 

Systems fl48p and (H^ are two different matrix generahzations of the Davey-Stewartson 
equation (DS-III-a and DS-III-b). In the scalar case equation ( l48l) was investigated 
The vector version of DS-III fHHj) and its Lax representation given by operators 
(jHJ and (H5|) in case Ci = 0, = 1 were introduced in^i"— . 

Let us consider (H7|) in the case where /3 = 1, it := q and M.q := ^ are scalars. Then 
(H7|) becomes 

a2^it2 = cuxx + ciUyy + VU + ficiSu, voy = -2^(|Mp)^, Sx = -2{\u\^)y. (50) 
Setting c = Ci = 1 and = 1, as a consequence of (150|) we obtain 

tt2^^t2 = ""xx + ^^j/s, + 'S'ltt, = -2(|up)^^ - 2(|wp)yy, (51) 

where Si = Vq + S. This is the well-known Davey-Stewartson system (DS-I) and (H7|) is 
therefore a matrix (noncommutative) generalization. The interest in noncommutative 
versions of DS systems (in particular, solution generating technique) has also arisen 
recently In Section [IV] we will consider dressing method for matrix generaliza- 

tion of DS system ( 147|) that leads to its exact solutions. In^^ DS-II system was also 
obtained from system ( H6l) and its Lax pair (Hlj), (l45l) after the change x ^ z, y ^ z. 

2. n = 2, ; = 2. 

M2,2 = M2,2 - C2(Lo)3 = a^dt, - 02/3393 - _ + 3/32c2q^7Wo^-'rT 
-3^C2qModyD-'r^qMoD-h^ + 3Pc2qMoD-'{r'^q}yMoD-'r^ (52) 
+3c2^2aj,qA^o^"^rTaj,. 
The commutator equation [Lq, M2,2] = is equivalent to the system: 

tt2qt2 - <lxx - C2(^''^qyyy - Voq + 3c2(3ciyMoSi + 3f3c2qMoSiy - 3c2qA^o'S'2 
-3c2qD-i {A^or^qA^o^i - MoSiMor^q} = 0, 

-a2r^ - rL - /3'c2rT , - r-^Vo + 3/3c25iAlorJ + 3c252-Mor^ (53) 

-3C2D-1 {^iMorTq - r^qAloSi} Mor'^ = 0, 

/3i;o, = -2(q>lor^)., ^i. = /3(rTq)„ ^2. = /3'(rjq),. 

13 



In this case the operator M3 i has the form: 

M3,i = M3,i - Ci(Lo)2 = asdt, - - v,D - vo - Ci/?^^^ + 2c,/3qMoD-'r^ + 
+2cil3(iMoD-^r^dy. 

(54) 

The equation [Lq, M3 1] = is equivalent to the system: 

«3qt3 = ^xxx + Ci(3'^qyy + Viq^^ + Voq + CiqA^o'S'i, 
-asrj.^ = -rj^^ + C2/3^rJy - (r'^fi)x + r^^^o + CiS'iA^or'^, 

Pvo,y = -3(q^A^or"^)x + [qMor'^ ,Vi], /3vi,y = -3(q7Wor^)x., S^.^ = -2/3(rTq)^. 



(55) 



n = 3, / = 2 



(57) 



-3/3c2qModyD-^r'^qMoD-^r^ + 3Pc2qMoD-^{r^q}yMoD-^T'^ (56) 
+3c2^2^j;q7Wo^"^r^5j,. 
The Lax equation [Lq, M3_2] = results in the system: 

«3qi3 - c<ixxx - C2f3^qyyy - cviq^ + 3c2/3qyXo'5'i + 3c2/3qA^o5'ij; - cvoq- 
-3c2qMoS2 - 3c2qD-^ {Mov'^qMoSi - MoSiMoT^q} = 0, 
a3rj - crj^^ - f3^C2ryyy - cr]^vi - cr'^vi^ + 3(3c2SiMor] + cr~^VQ+ 
+3c2S2Mor'^ - 3c2D-^ {^lA^or^q - r^qA^o^i} A^or^ = 0, 
I3vi^y = -3{qMor^)x, Si^ = f3{v~^q)y, 
(3voy = -3{q^Mor'^)x + [qA^or^, ^i], ^2^ = /3^{ryq)y. 

System (l57j) and its Lax pair Lq (jH]) and M3 2 (156|) was investigated in^. We list some 
reductions of this system: 

(a) a3,/?,c, C2 6 M, = q*, Aio = Ai^. The operators Li and M3 are then skew- 
Hermitian and (1571) takes the form 

a3qt3 - C^ixxx - C2f^^qyyy " CViq^ + 3C2/3qyMoSi + 3C2/3qMoSiy- 

-cvoq - 3c2qA^o'52 - 3c2qD"^ {Moq*qMoSi - MoSiMoq*q} = 0, 

(58) 

^Viy = -3(qAloq*)x, Si^ = (3{q*q)y, 

/3voy = -3(q^Aloq*)x + [qA^oq*,^^i], S2X = /3^(q^q)y 
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In the scalar case (A^ = m = 1), setting W 3 fj, := Aio, q{x,y,t3) := q{x,y,t3), 
(158|) reads 

"3% - cq^xx - CiQyyy + Sc/i^^ ^ \q\ldy + 3c2//gy j \q\ldx + 

+3c2/ig y {qqy)ydx + 3c^g y {qxq)xdy = 0. (59) 
In the real case q = q, (155]) becomes 

asqts - Cfixxx - C2(^''^q.yyy - CViq^ + SCif^C^yMoSi + 3C2qMoSiy- 

-cvsq - 3c2qMoS2 - 3c2qD"^ { A^oq^qA^O'S'i - A^O'S'iA^oq^q} =0, 

(60) 

I3viy = -3(qA^oq^)x, Six = /3(q^q)y, 

^Voy = -3(q^A^oq^)x + [qMoq^,Vi], 82^ = (3^{qyq)y. 
In the scalar case (A^ = m = 1), writing Aio = fi and q = q{x, y, t^) = q{x, y, t^), 
after setting y = x and c + C2 = — 1,/3 = 1 (I6OI) takes the form 



o^sQta + Qxxx - QfJ-q'^Qx = 0, (61) 

which is the mKdV equation. The systems (158|) and (l60l) are therefore, respec- 
tively, complex and real, spatially two-dimensional matrix generalizations of it. 

(b) /3 = 1, MqT^ = V with a constant matrix v. In terms of u := qz/, (1F7|) takes the 
form 



Q!3Mt3 - CUxxx - C2Uyyy + 30!) I MxC^yJ + 3C25j^ \^ "^^^ 

— c [u, fijfiy^ u — 3c2M [n, Sijfix^ = 0, 

^ j [u,vi]dy - 3c2 y [S'i,M](ix^ = 0, Viy = -3ux, Six = Uy. (62) 
In the scalar case (A^ = 1, m = 1), this reduces to 

astifg - CUxxx - C2Uyyy + 3cD I ^y Uxdy^ m| + 3c2dy ^'^ (^j Uydx^ } ^ ^^^^ 

which is the Nizhnik equation-^-. The system ( l62l) thus generalizes the latter to 
the matrix case. 

II. A; = 2 

Now we will consider two Lax pairs connected with the operator L2: 

L2 = P2dr, - - 2u - qMoD-^r'^ . (64) 
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1. n = 2, / = 1 



M2,i = M2,i - ci(L2)2 = a^dt, - - 2u - {(3ldl - 2(32dr,D^+ 

- 2(32Ur^ - AI32udr^ + AuD'^ + 4m^D + 2u^^ +2q^A^or^ + (65) 
+ 2ciMQY^D + - 2/32qA^o^"^5,2rT) . 

Lax equation [L2, M2,i] = is equivalent to the system: 

a2Ut^ - (32Ur2 = (qA^or^)x - Ci{L2{q}Mor^ + qMom{r}y)^, ^^^^ 
M2,i{q} = 0,Ml,{r} = 0. 

In vector case (A^ = 1) and under additional reductions Ci = 0, 0:2, /32 G zM, A^o = 
— A^Q, u = u and q = r system ( l66l) reduces to (2+l)-generalization of Yajima-Oikawa 
equation ( 12T]) . 



2. n = 3, / = 1. 



Ms.i = M3,i - ci(L2)2 = a3at3 - - 3tiD - I (m, + /32D-'{ur,} + qXor^) 
-ci {Pidl - 2^25.,D2 + /)4 _ 2/32 u,, - Ap2udr, + AuD^ + Au^D + 2m,, + 
+2q,A^orT + 2qMor^D + 4^^ _ 2P2qMoD'^dr,r^) . 

The equation [L2, Ma^i] = is equivalent to the following system: 
a?,Ut^ - - 3mm, - |/32(qAlor'^)r2 + 

+ |(qA^orJ - q,7Wor^)x + + ^2l^"H^^rJ + qA^or^]- 

-[qyWor^,^^] - CiL2{q}7Wor^ - CiqA^o(^5{r})^ = f/^I^^H^r^rJ, 
M3,i{q} = 0, M3^{r} = 0. 



(67) 



(68) 



In the vector case (A^ = 1) under additional reductions ^3 G M, /32 G i^i A^o = ~J^oi 
q = r and u = u system ( l68l) reduces to (2+l)-generalization of higher Yajima-Oikawa 
system (p3l) . 

Remark 4. After putting 02 = and 03 = m formulae ^^B^-^SBij we obtain examples of 
a new (l+l)-dimensional Matrix k-cKP hierarchy'^^ (see case 2 after the proof of Theorem 
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IV. DRESSING METHODS FOR EXTENSIONS OF (2+l)-DIMENSIONAL 
K-CONSTRAINED KP HIERARCHY 



In this section our aim is to consider hierarchy of equations given by the Lax pair fl30p . We 
suppose that the operators and M„ ^ in (I30p satisfy the commutator equation [L^, M„ = 
0. At first we recall some results from^. Let N x K-matrix functions ip and ip be solutions 
of linear problems: 

LkW} = ^A, Ll{^p} = A, A e MatK^Q- (69) 

Introduce binary Darboux transformation (BDT) in the following way: 

W = I -<f{C + D-\i;'^^}y^ D-^^'^, (70) 

where C is a i^T x i^T-constant nondegenerate matrix. The inverse operator W'-^ has the 
form: 

= I + ipD-^C + D-'{ij^^}y'ij^. (71) 
The following theorem is proven in^. 

Theorem 2.'^'' The operator := WLkW'^ obtained from Lk in (W\) via BDT ( fTOl) has 
the form 

k 
j=0 

where 

Mi = CA- A^C, $ = ipA-\ ^ = M"^'^, A = C + D-^itjj'^ip}, ^^^^ 
ci = W{q}, r = iy-i'"{r}. 
Uj are N x N -matrix coefficients depending on functions Lp, ip and Uj. In particular, 



,if{C + D-'{tP^ip}y'^^ . (74) 



Exact forms of all coefficients Uj are given in^. 
The following corollary follows from Theorem [2) 

Corollary 3. The functions $ = <pA-'^ = W{<p}C-^ and ^ = ^A'^'^ = W-^^^{iIj}C^'-^ 
satisfy the equations 

= $CAC-\ Ll{^} = ^C^AC^'-\ (75) 
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For further purposes we will need the following lemmas. 
Lemma 1. Let be a matrix of the form 

Mi+i = CA'+^ - (A^)'+^C, / e N. (76) 
The following formula holds: 

I 

Mi+i = J2 CA'C-^MiC-\A^y-'C. (77) 

s=0 

Lemma 2. The following formula 

I 

$A^i+iD-i^^ = J2 ^s]MiD-^'^^[l - s], (78) 

s=0 

holds, where 

m ■■= (L.y m ■.= {tiym. (79) 

Proof. Lemma 2 is a consequence of Corollary 1 and formula (177|) of Lemma 1. Namely, the 
following relations hold: 

s=Q s=0 

□ 

Now we assume that the functions ip and ip in addition to equations (169|) satisfy the 
equations: 

M^,i{^} = Q^A'+^ = ciLi+'{^}, M;,{^} = Q^A'+^ = ciiLiy^'W. (80) 
Problems flSU]) can be rewritten via the operator M„ ^ (15^ as: 

M^,iW} = o, m;,w = o. (81) 

The following theorem for the operators Mn^i (130|) and M„ / f l39|) holds: 

Theorem 3. Let N x K -matrix functions (p, ifj he solutions of problems ^U9\) and [ED^] . The 
transformed operator Mn,i := WMn^iW~^ obtained via BDT W ([T^j has the form: 

+ciY!s=o^[s]M,D-^^^[l -si A^ = Y.to^^D\ 
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where the matrix Ain o-nd the functions q, r, $[s], \E'[/ — s] are defined by formulae IT/^) . 
Ff^) and q[j], r[j] have the form 

m = (^i){q}, m = (^i)^{r}, (83) 

Vi are N x N -matrix coefficients that depend on the functions Lp, i/j and Vi . The transformed 
operator Mn,i = W Mn,iW~^ has the form: 

M^j = WM^jW-' = Mn,i - ci{Lk)'^\ (84) 

where is given by l[7^) . 

Proof. We shall rewrite the operator M„ ^ (130|) in the form 

n 

Mn,i = andt„ - J2 ^iD' - ciqMoD-'i^, (85) 

i=0 

where A^q is an m(/ + l) xm(/ + l)- block-diagonal matrix with entries of Aio at the diagonal; 
q := (q[0], q[l], • • • , q[^]), f := (i"[^], r[/ — 1], . . . , r[0]). Using Theorem [2] we obtain that 

n 

Mn,i = oindt^ - ^^D' - c^^o^^^'r^ + '^Mi+iD-'^^ , (86) 

i=0 

where q = W^{q}, q = VF"^'"^!?}. Using the exact form of q and f we have 

q = ^^{q} = (^^{q[0]}, • • • , W{<i[l]}), f = W-'^^{i} = {W''^^ {v[l]} , . . . , W'"'^ {v[0]}) . 

(87) 

We observe that 

W{q[i]} = WU{(i} = WVW-^{W{ci}} = L*{q} =: q[i]. (88) 

It can be shown analogously that iy~^''^{r[2]} = L'^'*{iy~^''^{r}} = L"^''!?} =: r[i]. Thus we 
have: 

I 

^iMoD-'i^ = J2^[j]MoD-H'^[l-j]. (89) 

j=0 

For the last item in fl86l) from Lemma 2 we have: 



^Mi+iD-^^^ = J2 ^s]MiD'^'^^[l - s]. (90) 

s=0 

Using formulae (186|) . ( 189|) . (190|) we obtain that the operator Mn^i has form (182|) . The exact 
form of the operator M„ ^ follows from formula ( 182|) and Theorem [2j □ 
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From Theorem [3] we obtain the following corollary. 

Corollary 4. Assume that functions and ip satisfy problems [6^) and / fgOj) . Then the 
functions $ = W{ip}C^^ and'^ = W~^''^{ifj}C~^'~^ (see formulae ITF^) ) satisfy the equations: 

M^,i{^} = M„,,{<|.} - ciiUY^'i^} = 0, m;,{vi/} = m;,{vi/} - QiLiY^'m = o, (9i) 

where the operators L^, Mn,i md M^^i are defined by H^l . ^SW '^'^^ W^ - 



As it was shown in Sections [ITllIIII the most interesting systems arise from the (2+1)- 
dimensional k-cKP hierarchy ([H]) and its extension (I30p after a Hermitian conjugation reduc- 
tion. Our aim is to show that under additional restrictions Binary Darboux Transformation 
W (170|) preserves this reduction. We shall point out that a differential dressing operator 
that was used in^^'^'^ does not satisfy such a property. It imposes nontrivial constraints on 
the dressing operator (see, for example,^-). The following proposition holds. 

Proposition 1. 1. Let ifj = and C = C* in the dressing operator W l[7I^) . Then the 
operator W is unitary (W* = W~^). 

2. Let the operator ^EW be Hermitian (skew-Hermitian) and Mn^i ^ES) be Hermitian 
(skew-Hermitian) . Then the operator = WL^W'^ (see ( 17^) ) transformed via the 
unitary operator W is Hermitian (skew-Hermitian) and Mn,i '■= W Mn,iW~^ is 
Hermitian (skew-Hermitian). 

3. Assume that the conditions of items 1 and 2 hold. Let A = A in the case of Hermitian 
Lk (A = —A in skew-Hermitian case). We shall also assume that the function ip 
satisfies the corresponding equations in formulae / fgPj) and ^8D\). Then Aii = —Ml^ 
(Ml = Ml) and = $ (see formulae ITI^) ). 

Proof. By using formulae ( 170]) and f lTT]) it is easy to check that W is unitary (in the case 
ip = Lp, C = C*). Assume that Lj. = LI (in the case of skew-Hermitian analogous 
considerations can be applied). Then we have (Lk)* = W^^'*LIW* = WLkW~^ = Lk- 
Analogously we can show that W maintains Hermitian (skew-Hermitian) property for Mn^i- 
The formulae from item 3 can be checked by direct calculations. □ 

Now we list several examples of dressing for Lax pairs and the corresponding equations 
from Section IIIII 
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Consider dressing methods for equations connected with the operator Lq Assume 
that and -0 are x X-matrix functions that satisfy the equations 



LoW} = ipA, Ll{iP} = tpK, Lo := (3dy, (92) 

By Theorem [2] we obtain that the dressed operator Lq via BDT W f lTD]) has the form 

Lo = WLoW~^ = Pdy + <l>MiD-^^^. (93) 

1. n = 2, / = 1. Assume that N x X-matrix functions ip and ip in addition to equations 
(|92|) also satisfy the equations 



y 



(96) 



M2W} = c.ipA^ = c,Ll{ip}, M^{iP} = ci^A^ = c,i^)^{^P}, M, := a2dt,-D\ (94) 

By Theorem [3] we obtain that the transformed operator M2 has the form 

M2 = VTMziy-^ = a2dt,-D^-vo + Lo{^}MiD-^^^ + ^MiD-\{Ll){-^}y. (95) 

By direct calculations it can be obtained that vq = 2{ipA-^tp'^)^, A = C + D'^ {iIj'^ ip} . 
It can be easily checked that 

= v?AA-V^ - /3v7A-iD-i{V'T^|^^-i^T _ (^A-^A^V)^ = 

= v?A-i(CA + D-i{V'^</'}A)A-V^ - /3v?A-iD-i{V'^V'}yA"V^+ 
+V9A-i(-A^C - A^D-^{ij~^p})A-^-i/j'^ = 
= pA-\CA + (5D'^{ij^py})A-^ij'^ - ^ipA-^D-^iij'^ p}yA-^ij'^ + 
+pA-\-A^C + l3D-^{'il;^p})A-^ij^ = <^Mi^^ . 

From the latter formula we obtain that 

^voy = 2/3((/;A-V^)., = 2{^Mi^^).. (97) 

From Corollary H] we see that the functions $ = LpA~^ and \E'= t/^A^'^^ where A = 
C + D~^{il:^ip} (see formulae (175]) ) satisfy equations (191 p . After the change q := $, 
r := ^, Mo ■= —Ml, Vo := -Oo from formulae f l^ and f p7|) we obtain that N x K- 
matrix functions q, r, an A x A^- matrix function Vq and a X x matrix Mq satisfy 
equations ( H6|) in the case c = 1. In the case of additional reductions in formulae ( p2|) - 
(jMj) : ^2 G zM, /3 G M, ci G M, A = -A, ip = p and C = C* in gauge transformation 
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operator W (170|) from Proposition [T] we obtain that the functions q := $ and vq = 
Vq = 2{ipA~^(fi*)x satisfy matrix DS system (H7j) in the case c = 1. 

Consider more precisely the exact solutions of a vector version {N = 1) of DS system 
(H7|) in the case c = 1. Assume that the K x X- matrix A in (p2|) is diagonal: A = 
diag(Ai, . . . , A/^), Aj G C. Let us fix arbitrary natural numbers Ki, . . ., Km- Denote 
by K the number K = Ki + . . . + K^- Our aim is to present the exact form of solution 
of m-component DS system (1471) (q = {qi, . . . qm)) Assume that a K x X- matrix C 
has the form: 



fJ-s 



±1, s = 1, m. 



(98) 



-^i + -^i/ i,j=i 

I.e., the matrix C has diagonal blocks Ck^ of dimension Kg on its diagonal, fig, 



s = 1, m, are numbers each equal to 1 or —1. It can be checked by direct calculations 
that the matrix M.i := CA + A*C has the form: 



Ml 



( HiEk^ Oki,K2 Oki,K3 
0k2,Ki fJ'2EK2 Ok2,K3 

Ok3,Ki Ok3,K2 IJ-sEks 



, 0K3,Km-l 



0k2,K„ 

0K3,Krr 



\ 



(99) 



OXm-uKi Ok^_^,K2 OXrn-uKs ■ ■ ■ fJ-m-lEK^n-l ^K^^^^Rrr, 

\ Ok^.Ki Ok^^K2 0Km,K3 ■ ■ ■ OKm,K^-i l^mEK^ J 

where by E^^ we denote the Kg x i^^-square matrix consisting of 1. I.e., E^^ = 
ij^lxs, where = (1, . . . , 1) is 1 x Xg-vector consisting of 1. Ok^^Kj is the 

Ki X Xj- matrix consisting of zeros. The matrix Aii can be rewritten as A^i = 
diag(/iil]^^lii-j, . . . , ^Umlx^lii'm)- It can be noticed that the matrix A4i admits the 
factorization: 



Ml = PaP*, P = (e^,, . . . ,eK™),cr = diag(/ii, . . . , [Xr. 



ex. 



(100) 



In formula ( llOOp by Ok. = (0, . . . , 0) we denote the 1 x Kj-vector consisting of zeros. 
ckj denotes the 1 x i^'- vector consisting of row vectors Okj and Ik^ 

Example 1. Consider the case m = 2, Ki = 1, K2 = 2. Then formulae ^9^-^0^) 
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for the matrix Ail become 



( 



Ml = 



\ 



/ii 

^2 ^2 

y iJ,2 J 



, Ml = P(tP* 



1 
y Oly 



//I 

/i2 



1 
1 1 



(101) 



We put /3 = 1 in (p2|) and choose solution of systems ( p2|) and (IMj) in the following 
form: 



ip = {(fii, . . . (Pk), Pj = exp <j ( ] t2 + TJjX + XjV J>, Aj, Vj G 



a2 



(102) 



Using Corollary m and formula W!\ we see that the functions $, = 2{ip/S. ^p>*) 
S = -2D-^{<^*<^}y = 2(A-i)y and Mi = PcrP* satisfy the matrix DS-system (gZl): 



a2^t2 = '^xx + ci%y + - Ci'^PaP*S, a2 G ^1 
voy = 2($PaP*<l>*),, ^, = -2(<1>*<I>),. 



(103) 



Define functions q, vq and S in the following way: 

q = ^P = pA-^P =p {C + D-'{ip*ip})-'p, z;o = 2(¥pA-V%., = 2P*(A-1),P, 

(104) 

where the matrices C and P are defined by (pS]) and fllOOp . The integral with 
respect to a: in the formula D~^{Lp*ip} (I104p is realized in the form 



Va + Uj / ■ ■ -, 
J l,J=l 



K 



(105) 



Then from equation fll03p it follows that functions fll04p satisfy the m-component DS 
system 

Oi2^.t2 = Qxx + Ciqyy + t;oq - CiqaS, a2 E iM, 
Voy = 2(qaq*)^, = -2{q*q)y, 

that can be rewritten as: 



(106) 



«2gi,t2 = qi,xx + ciQi^yy + voQi - ci Y^JLi f^jQjSji, i = i,m, 02 e n 

Voy = 2(EJLi/^iklj)^' Sji^^ = -2{qjqi)y, i,j = l,m. 



(107) 
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Remark 5. By the well-known formulae from the matrix theory solutions can 
be rewritten as 



det 



A CK, 

(p 



det A 



^0 = 2(¥pA-V): = -2 



/ /Av^*\\ 
det I 

V9 



det A 



=2(P.,)*(A-i),P., = -2 



det 



A P, 
fP.i)* 



det A 



V 



l,m; 



1, m. 



J 



(108) 



(109) 



(110) 



In formulae ^108\) - fTTW we set Kq := 0. By Sij, i,j = l,m, we denote the elements 
of the corresponding matrix-function S. By qj we denote the elements of the vector- 
function q = {qi, . . . , qm)- Pi denotes the i-th vector-column of the matrix P. 



In the case m = 1 system fll07p reduces to the following one: 



"25*2 = Qxx + Ciqyy + Voq - CifiiqS, 02 G iR, 
Voy = 2(/ii|g|2)^., = -2{\q\'^)y, 
and the formulae from Remark \5\ take the form: 



;iii) 



det 



q = 



I1 T 

K 



A ll 

If 



det A 



S = 2(lf,A-'l 



I1 T- 



^K)y 



Ik 
detA 



vo = 2(<^A-V*) 



det 



A if* \\ 


detA 
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and represent a i^-soliton solution of the scalar DS system. The latter in the case 
K = 1 takes the form: 

_ exp(g) ^ Re(Ai)exp(2Re(g)) _ Re(z/i) exp(2Re(^)) 
^ - ^ - ' " Re(Ai) ' ^^^^> 

where A = + ^^exp {2Re{e)) and 6 = (^%^) h + i^ix + \iy. We shall 

point out that DS equation fllll|) consists of two special cases: 

(a) /ii = 1. In this case formulae fllUSp -f lTTU]) from Remark |3] and flH2p represent 
regular solutions of DS-equation (lllip . 

(b) /ii = —1. In this case Remark |5] (formulae (I108p -( ITT0|) ) and formula ( I112p give 
us singular solutions of DS-equation fillip . 

In a similar way dressing technique has been elaborated for several multicomponent 
integrable systems from the k-cKP hierarchy 

n = 3, / = 2. Assume that in addition to equations fl92p the functions Lp and satisfy 
the equations: 

M^iv} = C2^A' = M^{tP} = C2^A^ = C2(L5)3{^}, M3 := a^dt, - D\ 

(113) 

As it was done in the case n = 2, / = 1, we obtained that the transformed operator 
M3 via BDT W ([70]) has the form 

2 

Ms = WM:,W-^ = a^dt,-D^ -v^D-v, + C2Y,L],mMoD-\{Llf-\^>}y , (114) 

where -Oi = 'i{LpA^'^il)^)xi Vq = SLpA'^ip'^ {(fA^^ifj^);j. — 3{ipx'^^^'^'^)x- After the change 
q := $, r := \E', Aio '■= —-Mi, Vi := Vi, vq := vq it can be shown that the functions q, 
r, Vi, Vq with the matrix A^q satisfy equations ( ISTP in case c = 1. 

Analogously, dressing methods can be done for the other Lax pairs from Section IIIII 
As an example we consider equation ( !68|) with Lax pair (167|) : 

A; = 2,n = 3,/ = 1. Assume that x K-matrix functions (p and if) satisfy the 
equations: 

L^W} = ^A, Llii;} = ^A, L2 := /^a^., - D', 

(115) 

MsW} = ciipK\ MJ{^} = ciiljk\ M3 := a^dt, - D\ 
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Using Theorem [2] we obtain that the dressed operator L2 has the form: 

L2 = p2dr^- D'^ -2u + ^MiD-^^^ , (116) 

where the functions $ and \1/ are defined by (I73|) : u = {(pA^^ilj^)^. Using Theorems 
[2] and |3l their corollaries and direct calculations it can be checked that the functions 
u = {ipA~^ifj'^)x, q := $, r = ^ and the matrix Mq := —Aii satisfy equations (168|) . 
We note that in the case Ci = 0, = 1 we obtain solutions for the (2+l)-dimensional 
generalization of higher Yajima-Oikawa equation: 

agMig - iw^^^ -3mm^ - 1/32 (qA^or'^)r2 + i (qA^or J -qx A^or'^)x=|/?|^"H«r2r2} , 
"sqtg - <ixxx - ^uq^ - l{ux + f32D-^{ur^} + qA^Qi"^) q = 0, 
a^Tt^ - r^xx - 3{ur)x + j {u^ + /32^~H^t2} + qA^or^) r = 0, 

(117) 

that was investigated as a member of (2+l)-dimensional k-cKP in^>^. Its solutions via 
differential gauge operators were considered in^*^. It should be mentioned that the 
solution generating technique via Binary Darboux Transformation W (170|) presented 
for matrix DS-system in item 1 can also be elaborated for system flllTp in the case of 
Hermitian conjugation reduction {u is a real-valued function, r = q, eta G M, /32 G i^, 
Ml = -Mo). 

We also point out that some of the matrix equations presented in items 1-3 and their solutions 
were also investigated recently—*^. 

V. CONCLUSIONS 

In this paper we proposed a new (2+l)-BDk-cKP hierarchy ( 130|) that generalizes (2+1)- 
dimensional KP hierarchy ([9]) which was introduced idM^ and rediscovered in^. For some 
members of (2+l)-dimensional k-cKP hierarchy (|9]) (e.g. ( 1211) and (123|) ). solutions were 
obtained via the binary Darboux transformation dressing method^ii^. Dressing methods 
for (2+l)-dimensional extension of the k-cKP and modified k-cKP hierarchy via differential 
transformations were elaborated recently irt22.New (2+l)-BDk-cKP hierarchy fl30p extends 
hierarchy ([9])^i^i^i^. It also contains integrable systems and their matrix generalizations 
that do not belong to (2+l)-dimensional k-cKP hierarchy (130|) . (2+l)-BDk-cKP hierarchy 
also extends matrix KP hierarchies such as (t^, r^) -hierarchy-^ and (7^1, (T5)-hierarchj4S. 
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Some members of new hierarchy (130!) such as Davey-Stewartson systems (DS-I,DS-II,DS-III), 
matrix generahzations of (2+l)-dimensional extensions of the mKdV and their representa- 
tions via integro-differential operators were considered recently in^. Dressing methods via 
binary Darboux transformations presented in Section HVl give us an opportunity to construct 
the exact solutions for equations that are contained in (2+l)-BDk-cKP hierarchy (1301) . The 
most interesting systems obtained from hierarchy (130|) arise after a Hermitian conjugation 
reduction. It is much more suitable to use BDT for dressing methods in this situation rather 
then a differential operator. Elaborated dressing technique method provides us also with a 
possibility to investigate exact solutions of matrix equations that (2+l)-BDk-cKP hierarchy 
(150]) contains. We shall note that the interest to noncommutative equations has also arisen 
recently^^"^'^ . Thus one of problems for further investigation is the generalization of (2+1)- 
BDk-cKP hierarchy to the case of noncommutative algebraa^i^"— (namely, the elements q, 
r, Uj and Vi belong to some noncommutative ring). It is also very important to point out that 
initial operators that we have chosen for dressing in ( l92l) . ( !94l) . ( 1113^ . ( lllSp are differential. 
However, the statements of the theorems in Section HVl (in particular Theorem [2] and [3]) ) con- 
cern initial integro-differential operators. The case of initial (undressed) integro-differential 
operators should also be elaborated. We recently obtained that for (l+l)-dimensional in- 
tegrable systems it leads to classes of solutions that do not tend to zero at infinities (e.g. 
finite density solutions). Generalizations of (2+l)-dimensional k-cKP hierarchy (130|) also 
provides us with problems for further investigations. In particular, using (2+l)-BDk-cKP 
hierarchy it becomes possible to investigate generalizations of modified (2+l)-extended k- 
cKP hierarchy (modified (2+l)-BDk-cKP hierarchy) and elaborate dressing methods for it. 
Some members of this hierarchy were recently considered in=-- . In particular, we investigated 
Lax integro-differential representations for the Nizhnik equation^^ and (2+l)-dimensional 
extension of the Chen-Lee-Liu equation^. Representations for some of those systems in 
the algebra of purely differential operators with matrix coefficients can be found in^. Using 
(2+l)-BDk-cKP hierarchy we also intend to generalize (2+l)-extended Harry-Dym hierarchy 
presented in^^. We recently obtained mult i- component equation with the following Lax pair 
contained in (2+l)-BDk-cmKP hierarchy (bidirectional generalization of (2+l)-dimensional 
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k-constrained modified KP hierarchy): 
Lo = dy- qMoD-^v^D, a G C, 

Ms = asdt, + cqD^ - C2dl - ScqViD^ - ScqVsD + 3c2qyMoD-^dyr^ D+ 

(118) 

+C2qMoD-^dyr'^dyD - 3c2qModyD-h^ DqMoD-^r^ D+ 
+3c2qMoD-^{v^ci,}yD-^Mor^ D, Co,C2 G M. 

where q and r are 1 x m- vector functions, u, Vi and ^3 are scalar functions, Aio is an 
m X m-constant matrix. If we impose additional reduction in flllSp : 

1 \ 

q = {aiQi, aiq2, tt2-D {u}, ^2), r = q, Mq = diag(/2, a), a = \ , ai, 02 G M 

-1 0) 

(119) 

where I2 is the identity matrix of dimension 2x2; qi, q2 and u are scalar functions. The 
system that is equivalent to the Lax equation [LqjMs] = generalizes (2+l)-extension of 
modified KdV equation that was investigated in^"^ (reduction C2 = 0, q;2 = in (I118p -( ITT9|) ) 
and the Nizhnik equation fl63l) (reduction ai = in flllQp ). 

The following remark shows connections between constrained KP hierarchies and Lax 
pairs with recursion operators. 

Remark 6. It should be also mentioned that the majority of Lax representations with recur- 
sion operators (the representations connected with bi-hamiltonian pairs of the corresponding 
integrable systems) involve integro- differential operators^^ . A part of such representations is 
contained in the constrained KP hierarchies (presented in Section III and IV) with additional 
(not only a Hermitian conjugation reduction) specific reductions. 

From the last remark it follows that one of the essential problems for further investigation 
is an adaptation of the dressing method (Theorem [2] and |3] from Section IV) for such Lax 
representations with the recursion operators. The work in this direction is in progress. 
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VII. APPENDIX. PROOF OF LEMMA 1. 

Proof. We will use the following recurrent formulae that can easily be checked by direct 
calculation: 

M2 = CAC-^Mi + MiC~A^C-\ (120) 

Mi+i = CAC-^Mi + MiC-^A^C - CAC-^Mi^iC-^A^C. (121) 
At first we will prove the following formula: 

k k 

Mi+i=J2 CA'C~^Mi-k+iC-\A'^)''-V-J2 CA'C-^Mi-kC-\A^)''-'+^C, k < 1-2, (122) 

s=Q s=l 

via induction by k. Assume that fll22p holds for some k < 1 — 2. We shall show that then 
this formula holds for A; + 1 < / — 1 using f ll2ip : 

Mi+i = Yl CA'C-^Mi-k+iC-^{A^f^'C - CA'C-^Mi-kC-^{A^f-'+^C = 

s=0 s=l 

= Y CA'+^C-^Mi-kC-\A^)''-'C + Yl CA'C-^Mi-kC-\A^f-'+^C- 

- Y CA'+^C-^Mi-k-iC-\A^f-'+^C - Y CA'C-^Mi-kC-\A^f+^-'C = 

s=0 s=l 

= Y CA'C-^Mi-kC-\A^f+^-'C + Mi-kC-\A^f+^C- 

i=l 

k+1 _ fe+1 

- Y CA'C-^Mi-k-iC-^{A^f''^^C = Y CA'C-^Mi-kC-\A^)''-'+^C- 

i=l i=0 

- Y CA'C-^Mi-k-iC-\A^)''-'+\ 

i=l 

Thus, we have proven (I122p . After the substitution of A; = / — 2 in (177|) and using (I120p we 
obtain: 

Mi+i = Y CA'C-^M2C-\A^y-^-' - Y CA'C-^MiC-\A'^y-' = 

s=0 s=l 

= Y CA'+^C-^MiC-\A'^y-^-' + Y CA'C-^MiC-\A^y-'C- (124) 
- Y CA'C'^MiC-\A'^y-'C = Y CA'C'^MiC-\A^y-'C. 

s=l s=0 

This finishes the proof of formula ( !77|) and Lemma 1. 

□ 
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